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Abstract 

Spacetime supersymmetry is widely believed to play an important role in 
most fundamental theories of physics, and is usually invoked in order to address 
problems of naturalness. In this paper, we examine the question of whether 
supersymmetry itself is "natural" (i.e., likely to exist as a fundamental compo- 
nent of nature at high energy scales). Our approach to answering this question 
is based on a statistical examination of the heterotic string landscape, and 
our conclusion is that supersymmetry is an exceedingly rare phenomenon. We 
also find that the likelihood of supersymmetry appearing at the string scale is 
dependent on the gauge symmetries present at the string scale, with certain 
gauge groups strongly favoring the appearance of M=l supersymmetry and 
others not. This article summarizes several recent papers, yet also contains 
some new results. In particular, one new result is that the heterotic land- 
scape appears to favor either the non-supersymmetric Standard Model or an 
M=l SUSY GUT gauge group at the string scale; by contrast, the opposite 
outcomes (namely the MSSM or a non-supersymmetric GUT) are significantly 
disfavored. 
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1 Introduction 



Most theoretical frameworks for physics beyond the Standard Model involve the 
introduction of supersymmetry (SUSY), and there are many reasons why this is 
so. First, supersymmetry solves the technical gauge hierarchy problem. Second, 
supersymmetry provides a dynamical mechanism for triggering electroweak symmetry 
breaking. Third, supersymmetry improves the accuracy of gauge coupling unification, 
and fourth, it provides a dark matter candidate. As a result, supersymmetry is truly 
ubiquitous in particle physics, appearing virtually everywhere — except the data. 

However, through the years, lots of competing or alternative theories have been 
proposed. Some involve large extra dimensions and some involve small extra di- 
mensions. Others do not explicitly involve extra dimensions at all, yet contain new, 
strongly coupled sectors. Likewise, especially over the past decade, the substance of 
phenomenological model-building has changed dramatically. Indeed, it is now quite 
common that a talk introducing a new particle-physics scenario will begin with a 
litany of sequential assumptions that would have sounded increasingly fantastic to 
the ear of a physicist a mere decade ago. We are made of open strings. And we 
live on a brane. And the brane lives in extra dimensions. And the brane is wrapped 
and intersects other branes. And the extra dimensions are warped. And the warp- 
ing is severe and forms a throat. And the brane is falling into a throat. And so 
forth, and so on. Indeed, such scenarios now often form the backbone of cutting-edge 
model-building. 

Admittedly, all of this may sound highly unnatural, and it is excusable to yearn 
for the simpler days of the MSSM and their cousins, the SUSY GUTs. But is SUSY 
itself truly natural? What does it mean to be "natural" , anyway? 

There are many different notions of naturalness that have appeared in the litera- 
ture. For example, Dirac naturalness stipulates that an effective field theory (EFT) 
is natural if the dimensionless coefficients for all operators are ~ Oil) - - i-e., no 
exceedingly small or large numbers are allowed. In this sense, the large electroweak 
gauge hierarchy is unnatural, which is one of the biggest motivations for supersym- 
metry. Another notion is 't Hooft naturalness: even if such a number is small, it can 
be viewed as "natural" if its smallness is protected by a nearly unbroken symmetry. 

But neither of these addresses the question of whether a theory, even if "natural" 
in the above sense, is likely to be right. How likely is SUSY to be the correct theory? 

The word likely often causes us to shudder. Indeed, even though we constantly 
judge theories this way, we don't say this word aloud because the question of theoret- 
ical likelihoods seems more philosophical than scientific, especially when we have no 
data upon which to base our assessments. How likely relative to what? To all other 
theories that one can imagine? And who is doing the imagining? One might get very 
different responses depending on the identity of the unlucky proponent. Ultimately, 
we seem to be faced with a dead-end question. How can one compare the likelihood 
of one theory against another? 
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String theory provides a framework in which this question can be addressed in 
a mathematical way. This is because string theory provides a large set of possible 
vacuum solutions ("vacua", or string "models", collectively called "the landscape"), 
each of which corresponds to a different alternative universe with different physical 
laws. In this context, we can then place our likelihood question on firmer statistical 
footing, as originally advocated in Ref . pQ . In the landscape of possible string solu- 
tions, how many of these solutions are super ■symmetric? Is SUSY "natural" on this 
landscape, or relatively rare? 

This is the subject of this article. 

2 Our Study: the Heterotic Landscape 

In order to approach this question in a somewhat tractable way, we shall focus on 
the landscape associated with perturbative four- dimensional heterotic strings which 
are realizable using free-field constructions (such as those based on free worldsheet 
bosons or fermions). Our investigation will then focus on the fundamental question 
of determining the probability that such theories exhibit spacetime supersymmetry 
at the string scale. Note that this is different from previous discussions [2] in which 
a certain intrinsically supersymmetric framework [3] for the high-scale theory is as- 
sumed, and in which one asks about the likelihood that dynamical supersymmetry 
breaking is subsequently generated near the electroweak scale. 

In a recent paper [I], we conducted a random exploration of the landscape asso- 
ciated with such models, using the free-fermionic construction method [5J. In this 
construction, each model can be described in terms of its left- and right-moving 
worldsheet conformal field theories. These conformal field theories consist of ten- 
sor products of non-interacting, free, complex fermionic fields. Different models are 
then achieved in this method by varying the boundary conditions of these worldsheet 
fermionic fields around the two non-contractible loops of the worldsheet torus. The 
sets of allowed boundary conditions for the fermionic fields are restricted by numerous 
string self-consistency conditions which also must be applied. These self-consistency 
conditions guarantee that the string partition function can be viewed as the trace 
over a Fock space corresponding to a self-consistent string model. Other perturba- 
tive self-consistency constraints (such as conformal and modular invariance) are also 
imposed. 

This study was similar to earlier studies performed in Refs. [6j [7], and resulted 
in a data set of approximately 10 7 distinct self-consistent four- dimensional heterotic 
string models, each of which is tachyon-free and hence stable at tree level. To the 
best of our knowledge, this is the largest set of distinct heterotic string models ever 
constructed. However, this study has some important limitations. For example, this 
study examined string models for which all real worldsheet fermions could be grouped 
together to form complex worldsheet fermions. This effectively reduces the number of 
consistent sets of fermion worldsheet boundary conditions, but allows for the utiliza- 
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tion of many time-saving computational algorithms. The major phenomenological 
consequence of this restriction is that these string models have gauge groups with a 
fixed rank (namely twenty-two in four dimensions). 

In our study, two models were considered distinct if their spacetime phenomenolo- 
gies differ in some way. For example, two models were considered distinct if they differ 
in their amounts of unbroken supersymmetry, or in their gauge group or massless par- 
ticle spectrum. Moreover, using data gathered from our sample set of models, a novel 
technique [8] was utilized in order to extract stable statistical results which do not 
vary as a function of sample size. As discussed in Ref. [8], this method transcends a 
mere statistical analysis of the models in our limited sample set, and yields results 
which are likely to be indicative of the corresponding landscape associated with such 
models as a whole. The techniques by which such results are extracted is discussed 
in Ref. [8]. 

3 Results 

In this section, we shall describe the results of our study. As indicated above, 
the statistical methodology utilized in extracting these results is discussed in Ref. [8], 
and a full discussion of the results of Sects. 3.1 and 3.2 can be found in Ref. [3]. 

3.1 Prevalence of SUSY on the heterotic landscape 

One of the first questions addressed in our study concerned the prevalence of 
spacetime supersymmetry in the heterotic landscape. In other words, what percent- 
age of the heterotic string landscape consists of string models exhibiting an unbroken 
spacetime supersymmetry at the string scale? Clearly, answering this question is the 
first step towards addressing the overall issue of the naturalness of supersymmetry. 

Our results [3] are reproduced in Table [TJ Note that in quoting these results, we 
are explicitly focusing on only that portion of the heterotic string landscape which 
is stable at tree level. In other words, we are explicitly disregarding those non- 
supersymmetric portions of the landscape (which amount to approximately 32.1% of 
the total) for which tachyonic states exist at tree level. 

Tabled] represents our final partitioning of the tree-level four- dimensional heterotic 
landscape according to its degree of supersymmetry. There are several rather striking 
facts which are evident from these results. 

• First, we see that more than half of the stable heterotic landscape is non- 
supersymmetric and yet tachyon-free. Indeed, this proportion remains near half 
of the total even when the non-supersymmetric tachyonic models are included. 

• Second, we see that the supersymmetric portion of the heterotic landscape ap- 
pears to account for less than one-third of the full tachyon-free four-dimensional 
heterotic landscape. 
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SUSY class 


% of hptprotic landscaoe 


Af—O (tachyon-free) 


68.48 


Af=l 


30.78 


Af=2 


0.74 


A/"=4 


0.0044 



Table 1: Classification of the four-dimensional heterotic landscape as function of the number 
of unbroken spacetime supersymmetries. We are explicitly focusing on string models which 
are stable (and thus tachyon-free) at tree level. 

• Finally, models exhibiting extended (Af > 2) supersymmetries are exceedingly 
rare, together representing less than one percent of the full landscape. 

Of course, we stress once again that these results hold only for the tree-level 
landscape, i.e., models which are stable at tree level only. A priori, it is not clear 
whether these results would persist after full moduli stabilization. However, it seems 
likely that they would, given that most modern methods of moduli stabilization 
(fluxes, superpotentials, etc.) tend to further break (rather than restore) spacetime 
supersymmetry. Indeed, under these assumptions, these results then lead to a number 
of interesting conclusions. 

The first conclusion is that the properties of the tachyon-free heterotic landscape 
as a whole are statistically dominated by the properties of string models which do 
not have spacetime supersymmetry. Indeed, the A/"=0 string models account for 
nearly three-quarters of this portion of the heterotic string landscape. The fact that 
the J\f=0 string models dominate the tachyon-free portion of the landscape suggests 
that breaking supersymmetry without introducing tachyons is actually favored over 
preserving supersymmetry for this portion of the landscape. Indeed, we expect this 
result to hold even after full moduli stabilization (as has been argued within the 
context of Type I strings j9j), unless an unbroken supersymmetry is somehow restored 
by stabilization. 

The second conclusion which can be drawn from these results is that the super- 
symmetric portion of the landscape is almost completely comprised of A/"=l string 
models. Indeed, only 2.4% of the supersymmetric portion of the heterotic landscape 
has more than Af=l supersymmetry. This suggests that the correlations present for 
the supersymmetric portion of the landscape can be interpreted as the statistical 
correlations within the Af—1 string models, with the Af=2 correlations represent- 
ing a correction at the level of 2% and the J\f=4 correlations representing a nearly 
negligible correction. 

In fact, the SUSY faction of the full string landscape may be even smaller than 
quoted here. One reason is that free-field string constructions (such as we are employ- 
ing here) probably tend to artificially favor models with unbroken supersymmetry. 
Second, even when stabilized string models nevertheless exhibit spacetime SUSY at 
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the string scale, there remains the difficult quesion of determining the statistical like- 
lihood that this SUSY will survive all the way down to the electroweak scale [2J. At 
the very least, this should restrict the number of string models leading to weak-scale 
SUSY still further. 

Thus, we conclude that weak-scale SUSY is rather unnatural from a string land- 
scape perspective. On the one hand, this result shifts the burden of proof onto the 
SUSY enthusiasts, which represents a dramatic reframing of the underlying question 
of whether SUSY should exist at or above electroweak scale. But the conclusion 
that weak-scale SUSY is unnatural should not necessarily be viewed as a problem 
for string phenomenology. In fact, this result might even be considered good news: 
it implies that we will actually learn something about string theory and its pre- 
ferred compactifications if /when weak-scale supersymmetry is actually discovered in 
upcoming collider experiments! 

3.2 Correlations between SUSY and gauge groups 

We shall now examine the effects of supersymmetry on the probability of realizing 
different gauge group factors in this landscape. What percentage of string models 
with a given level of supersymmetry will contain a given gauge group factor amongst 
its unbroken gauge group at the string scale? We shall also be interested in knowing 
the likelihood of realizing different gauge group factors for the full landscape as a 
whole. 

Our results [I] are presented in Table [2j As can clearly be seen, supersymmetry 
has a profound effect upon the prevalence of different gauge group factors. Moreover, 
even independently of SUSY, there are some general trends which emerge from these 
results. These trends include: 

• A preference for SU(n + 1) over SO(2n) groups for each rank n. Even though 
these two groups have the same rank, it seems that SU groups are more common 
than the SO groups for all levels of supersymmetry. 

• Groups with smaller rank are much more common than groups with larger rank. 
Once again, this also appears to hold for all levels of supersymmetry. 

• Finally, the gauge-group factors comprising Standard-Model gauge group 
G SM = SU(3) x SU{2) x £7(1) are particularly common, much more so than 
those of any of its grand-unified extensions. 

3.3 SUSY naturalness 

At this point, we have presented both the unrestricted probability of finding differ- 
ent levels of supersymmetry on the heterotic landscape, and the restricted probability 
of finding different gauge group factors when a certain degree of supersymmetry is 
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Table 2: The percentage of heterotic string models exhibiting specific gauge group factors 
as functions of their spacetime supersymmetry. Here SU>§ and SO > \q collectively indicate 
gauge groups SU(n) and SO(2n) for any n > 5, while J\f refers to the number of unbroken 
super symmetries at the string scale. Note that the 7V=0 models are all tachyon-free. The 
rightmost column of this table is derived from the other columns using the landscape 
weightings in Table [H 



assumed. However, a more useful quantity might be the "inverse" of this last prob- 
ability, namely the probability of finding different levels of supersymmetry given a 
specific gauge group factor. This would give an indication of how "natural" each 
level of supersymmetry is for each different possible gauge group factor. 

In order to derive these probabilities, we can utilize the results presented above. 
First, let us recall some basic elements of probability theory. The results in the 
first four columns of Table [2] are necessarily conditional probabilities: each number 
in these columns represents the probability of finding a specific gauge group factor 
given a certain level of supersymmetry. Thus, if A represents the occurrence of a 
specific gauge group factor and B represents the occurrence of a specific level of 
supersymmetry, then the results presented in Table [2] are all of the form 

P(A I B) s f£™> . (3.1) 

Of course, what we now seek is not P(A \ B), but the "inverse" P(B \ A). In general, 
the relationship between P(B | A) and P{A \ B) is given as 

P(B I A) = P(A I B) . (3.2) 

Fortunately, all of the probabilities needed within the right side of Eq. (13. 2p are 
present in Tables [I] and [2j Specifically, P(A) is given within the rightmost column of 
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Table 3: The percentage of string models with different levels of supersymmetry as a 
function of different gauge group factors. Thus, if we know that a given string model gives 
rise to a specific gauge group factor at the string scale, this table lists the corresponding 
probabilities that this model will have various levels of unbroken supersymmetry. This 
table can therefore be viewed as the "inverse" of Table El 

Tabled while P(B) is given in Tabled] and P(A \ B) is given in the rest of Table [21 
Thus, given the presence of specific gauge group factor at the string scale, we can 
now determine the corresponding probability of finding different levels of unbroken 
supersymmetry at the string scale. Our results are presented in Table [31 Note that 
in calculating the results in Table [31 we have retained more significant digits than 
are explicitly shown in Tables [T] or [21 

By comparing these probabilities to the probabilities given in Table [U it is pos- 
sible to determine which gauge group factors tend to favor different levels of super- 
symmetry beyond their expected representations on the landscape as a whole. For 
example, we see from Table [1] that 68.48% of the stable heterotic landscape is non- 
supersymmetric. Thus, if a given gauge group factor is associated with models of 
which fewer than 68.48% are non-supersymmetric, then this gauge group factor can 
be said to preferentially favor supersymmetry. 

These results have a number of dramatic implications. 

• First, we observe that gauge group factors with large rank (greater than four) 
actually favor the appearance of unbroken supersymmetry. 

• Second, we observe that the gauge group factors which comprise the Standard 
Model gauge group do not generally favor supersymmetry. 

• Finally, we see that the SU{n) gauge group factors (with n > 5) and the excep- 
tional gauge groups £6,7,8 overwhelmingly favor Af = 1 supersymmetry. This 
preference is substantial, resulting from the combined effects of the individual 
probabilities contributing to Eq. (13. 2p . 

We thus conclude that the heterotic string landscape appears to favor either the 
non-supersymmetric Standard Model gauge group or an Af = 1 SUSY GUT gauge 
group at the string scale. However, the opposite outcomes (namely the MSSM or a 
non-SUSY GUT gauge group) are significantly disfavored. 
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One important caveat is that the gauge group factors presented in Table [3] do not 
generally specify the gauge group fully. Indeed, these gauge group factors could be 
part of either a hidden sector or the visible sector. However, the gauge group factors 
listed in these tables are necessarily among those which are explicitly present in the 
full gauge group of the string model at the string scale. 

4 Discussion 

In this paper, we have presented results concerning the prevalence of spacetime 
supersymmetry at the string scale and its possible statistical correlations with the 
unbroken gauge group which might also appear at the string scale. Since these 
two quantities (spacetime supersymmetry and unbroken gauge group) are completely 
independent in an ordinary quantum field theory based on point particles, these sorts 
of correlations represent true predictions of string theory and thereby provide one 
possible route towards answering the question as to whether supersymmetry is truly 
"natural" as a component of Beyond-the-Standard- Model physics. 

As we have seen, the results of this calculation show that spacetime supersym- 
metry is not generically a feature of the low-energy limit of string theory. However, 
spacetime supersymmetry is actually statistically favored for certain gauge group 
factors. 

There are some inherent limitations to these results which must continually be 
borne in mind. For example, these string models are generally unstable: the non- 
supersymmetric models generically have non-zero dilaton tadpoles, and the super- 
symmetric string models have flat directions. Thus, one might think that these 
results might change if only fully-stabilized models are considered. Unfortunately, 
no fully stable perturbative heterotic string models have ever been constructed. One 
could even argue that if the non-supersymmetry string models were required to be 
as "stable" as the supersymmetric string models (e.g., only have some finite number 
of flat directions), then these results would also change. However, at this point in 
time, the string models considered in this study are state-of-the-art and are as stable 
as many of the other classes of string models which have been considered in other 
statistical studies. 

Another issue facing this study concerns the extent to which these sorts of statis- 
tical correlations can be trusted, given that the full heterotic landscape has not been 
surveyed. However, this issue has been discussed in Ref. [8], and the methods devel- 
oped there have been used in order to extract each of the results quoted here. Thus, 
to the best of our knowledge, the results quoted here are independent of the size of 
our sample of heterotic string models, and thus should persist across the landscape 
as a whole. 

Finally, one could argue that the construction method utilized in this study neces- 
sarily only probes certain sections of the heterotic landscape. While this is true, these 
sections are the ones most likely to contain string models realizing non-abelian gauge 
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symmetries. As such, these are the sections most likely to give rise to non-trivial low 
energy phenomenologies. 

The interpretation of these results is also open to some debate. The probability- 
based definition of naturalness used in this paper is not the traditional one, and may 
only hold relevance for a landscape study such as the one we are performing. However, 
this definition of naturalness has the advantage of being applicable in a wide variety 
of contexts, and does not resort to any aesthetic or theoretical prejudices concerning 
the parameters that appear in effective Lagrangians. As such, probability-based 
definitions of naturalness may have inherent advantages over other definitions. 

There are several extensions to these results which are currently under investiga- 
tion. For example, we would like to understand how the presence of supersymmetry 
affects the statistical appearance of the entire composite Standard- Model gauge group 
Gsm = SU(3) x SU(2) x U(l), and not merely the appearance of its individual fac- 
tors. We would also like to understand how the presence or absence of supersymmetry 
affects other features which are equally important for the overall architecture of the 
Standard Model: these include the appearance of three chiral generations of quarks 
and leptons, along with a potentially correct set of gauge couplings and Yukawa 
couplings. This work will be reported elsewhere [10J. 
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